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Deseription 

Field ot the Invention 

5 The present invention relates to methods of and apparatus for, determining the location of an object and in partic- 

ular, but not exclusively to methods and apparatus which employ a magnetic field which is sensed at the object. 

Background of the Invention 

10 It has been long appreciated that If the magnetic field around a field generating element, for example a generatirig 
coil, can be accurately mapped then It might be possible to determine the location of a field sensor, for example a 
sensing coil, relative to the generating coil, from the signals sensed by such a sensing coil. However, a problern as- 
sociated with doing this is that there are in general many locations and/or orientations of the sensing coll within thp 
field of the generating coil that will provide the same characteristic sensing signals in the sensing coil. In order to use 

IS a magnetic field for this purpose, additional infonmation must therefore be provided. 

Prior art approaches to providing the additional information required comprise either moving the generating and 
sensing coils relative to each other, or scanning the axis of the generated field past the sensing coll. 

An example of the first approach is taught in US 3,644,825 wherein a system is disclosed for locating the position 
of a field sensor, comprising two orthogonal sensing coils, relative to a field generating element which relies on having 

20 knowledge of the direction of motbn of the sensor relative to the generator. It should be noted that this system cannot 
detect the location of an object unless there is such relative nrK>tk3n. and its direction is known. 

The second approach of scanning the axis of the generated field is disclosed, for position location In two dimen- 
sions, in US 3,121,228 and for position location in three dimensions in US 3,868,565. 

US 3,121,228 describes how the distance and direction of a sensor, again comprising two orthogonal sensing 

2S coils, relative to a fiek) generator, also comprising two orthogonal coils, can be determined. The two orthogonal gen- 
erating coils are driven in phase quadrature so that the axis of the resultant field Is caused to rotate within a plane. If 
' the sensor Is located within this plane then the axis of the fiekl is guaranteed to scan past the sensor, and. because 
at any given distance from a field generator the field strength will be a maximum at the field axis, the sensor will detect 
a maximum in field strength at this time. The voltage Induced in any one of the two coils forming the sensor will be 

30 dependent on the orientatbn of the coil relative to the field generator, and it is for this reason that in '228 two orthogonal 
coils are utilised In the sensor. The sum of these two voltages gives an indication of the distance between the sensor 
and generator, while the phase difference between the two voltages gives an Indication of the direction of the generator 
relative to the sensor. It is thus essential to the operation of the locatbn system of *228 that the axis of the field rotates 
and that two coils are present in the sensor. 

3S In US 3, 868, 555 this approach of scanning the axis, or maximum intensity vector, of the field past the sensor Is 

extended to allow locatbn of the sensor in three dimensbns. Whereas in two dimensions it Is sufficient merely to rotate 
the axis of the field within the plane to be sensed to guarantee It passing through the sensor. In three dimensions the 
axis would have to be rotated so that it described the surface of a sphere in order to be certain it encountered the 
sensor. To ensure that the axis passed through all points on the surface of a sphere the motion of the axis would be 

40 such that it encountered the sensor only very infrequently, and thus measurements by the sensor of the maximum field 
strength would also be infrequent. To avoid this the location system of '565 drives the generator coils in a complex 
fashion so that the field axis tracks and rotates around the position of the sensor. 

In order to locate the position of the sensor In three dimensions, according to the method of '565, three mutually 
orthogonal generating coils and three mutually orthogonal sensing coils are required and the three generating coils 

45 must be driven simultaneously by the three drive currents having amplitude and phase relationships between them 
which are controlled so as to direct the field axis towards the sensor. 

The approach taken in '565 further requires that the various equations governing the voltage induced in a sensing 
coil located and orientated in a particular alternating magnetic field are solved dynamically in real time I.e. during the 
acquisition of data from the sensing coil. This requirement, in additton to limiting the speed at which the sensor can 

so move while still being located successfully by the system, also means that should It be desired to locate more than 
one sensor, all apparatus will need to be duplicated for each additional sensor. 

US 4,710,708 discloses a positbn locatbn system, in which it is not necessary to scan the field axis. 708 employs 
multiple coil field generators and a single coil sensor, but utilises standard iterative algorithms to solve for all the var- 
iables of the relevant simultaneous equations, in a computationally intensive manner. 

55 

Summary of the Invention 

According to a first aspect of the present invention there is provided a method of determining the location and/or 
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orientation of a field sensor relative to a plurality of field generators of known location, each field generator comprising 
a plurality of collocated field generating elements, each generating element being energized and a measurement being 
made at the field sensor of the respective Held generated by the generating element, the method being characterised 
by the steps of: 

5 

1) for each field generator calculating, from the measurements of the field generated by each of its generating 
elements and an estimate of the orientation of the sensor, an estimate of the distance from that particular field 
generator to the sensor, and 

2) utilising the estimates of distance from step 1), and the known location of the field generators to calculate an 
10 estimate of the location of the sensor; 

3) employing the estimated location of the sensor from step 2) and the measurements of the field at the sensor to 
calculate a new estimate of the orientation of the sensor; 

4) repeating steps 1) to 3) iteratively. with step 1) employing the new estimate of sensor orientation from the 
preceding step 3), to improve the estimates of location and orientation of the sensor. 

15 

The method of the first aspect of the present invention thus enables the location and orientation of a sensor to be 
determined without either relative motion between the sensor and the field generating element, or scanning of the axis 
of the field. 

Furthermore, by calculating an estimate of the distance of the sensor from each field generator, a surprisingly 
20 accurate estimate of the position of the sensor Is achieved in a computationally simple manner. 

Since the method dissociates the stages of acquisition of data from the sensor, and processing of that data, rapid 
determination of the sensor location is facilitated. Furthermore the location of additional sensors may be detemnined 
simply by simultaneous measuring the field, generated by each generating element, at these other sensors and inde- 
pendently calculating their distances from the field generators. It should be noted that no modification of the field 
25 generating apparatus or method of driving the apparatus Is required In order to determine the location of a plurality of 
sensors. 

The applicants have discovered that advantageously the method of the first aspect of the present invention also 
allows the location of a sensor comprising a single sensing element, for example a sensing coil, to be determined, as 
will be explained subsequently. This is particularly advantageous for posrtioning applications in which two or more 
30 mutually orthogonal sensing coils, as required by prbr art techniques, cannot be used 

According to a second aspect of the present invention there is provided a method of determining the location of a 
field sensor, comprising a plurality of collocated field sensing elements, relative to a field generator, comprising a 
plurality of collocated field generating elements, each field generating element being energized in turn and the value 
of .the resultant field strength being measured at the field sensor, the measured field strength being dependent on the 
$s locatfon of the sensor within the respective field, the method being characterised by the steps of: 

1) calculating, by utilising all the values of measured field strength and an estimate of the direction of the sensor 
from the field generator, a direction dependent weighting factor for each field generating element so that the cal- 
culated field strength B is equal to the field strength B that would exist at the sensor if the axis of the field were 

40 directed towards the sensor, 

2) iteratively altering the direction dependent weighting factors to maximise B and thus to determine to a desired 
level of accuracy the direction of the sensor frorn the field generator, and 

3) employing the measured values of the field strength to calculate the distance of the sensor from the field gen- 
erator and hence, from the directbn of the sensor in step 2). the location of the sensor relative to the field generator. 

45 

This aspect of the invention thus provides a method of locating a sensor relative to a single field generator. 
The invention further provides apparatus suitable for carrying out the methods of the first two aspects' of the in- 
vention. 

Preferred embodiments of the present inventton will now be described, by way of example only, with reference to 
so the accompanying figures in which: 

Figure 1 shows a first embodiment of the invention, 

Figure 2 shows the cartesian coordinate system employed for a sensor of arbitrary orientation located at point P, 
Figure 3 shows schematically the resolution of the magnetic flux density at a sensor, 
55 Figure 4 shows the coordinate system employed to locate a sensor relative to a field generator. 

Figure 5 shows schematically a simulated circle of constant induced voltage, in a sensor, in two dimensions, which 
is employed in the first embodiment of the inventbn. 

Figure 6 shows schematically three simulated spheres of constant induced voltage each centred on a field gen- 
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erator, which is employed in the first embodiment of the invention, 

Figure 7 shows a flow chart of a first positioning algorithm used in the first embodiment of the invention. 
Figure 8 shows the same schematic as Figure 6 when the Iccatlon and orientation of the sensor have been deter- 
mined, 

5 Figures 9. 10 and 11 schematically show a coordinate system emptoyed in a second positioning algorithm used 

in the first embodiment of the inventton, 

Figure 12 shows images of an endoscope obtained using the positioning system of the present invention on the 
left, and inoages obtahed conventbnally using X-rays on the right, (a) shows a sigmoid loop, (b) shows an alpha 
loop, and (c) shows a reverse alpha loop, 
10 Figure 13 shows images of an endoscope within a patient obtained using the present positioning system on tKe 
left, and obtained using conventional X-ray imaging on the right, (a) shows an anterior view, and (b) shows a lateral 
view. 

Figure 1 4 shows a second embodiment of the Inventbn. 

IS. In a first embodiment the invention enables a sensor comprising a single sensing coil to be located in three dimen- 
sions relative to a plane defined by three field generators. 

Detailed Descrtptlon of the Drawings 

20 With reference to Figure 1 , three field generators 1 are mounted at known locations on a plane surface 2, Each 
field generator 1 comprises three electrically separate coils of wire (generating coils) 3 wound about a cuboid wooden 
former 4, which is approximately 40mm along one side. The three coils of each field generator are wound so that the 
axes of the coils are mutually perpendicular. The nine generating coils are separately electrically connected to an 
amplifier 5 which is able, under the direction of a controller 6, to drive each of the coils individually Each coil comprises 

25 40 tums of 0.45mm copper wire and has an inductance of approximately 75mH. 

The sensor 7 comprises a single sensing coil of 200 tums of 42swg wire on a ferrite core of diameter 0.6mm, and 
length 12mm. Larger sensor coils will in general be more sensitive to the electro-magrtetic fields generated by the 
generating coils, however the size of the coil is normally governed by the particular position location problem which is 
being addressed and frequently small sensor coils will be required. For an air-cored coil the sensitivity of the sensor 

30 depends on the area of the coil, however the sensitivity can be increased by utilising a high magnetic permeability 
material in the core, and in this case the sensitivity will depend more strongly on the length of the coil than on its 
diameter. The sensing coil is electrically connected to a measurement unit 8 which in turn is connected to the controller 
6. The measurement unit 8 comprises an analogue to digital converter, and a matched filter (not shown). 

In use, the controller 6 directs the amplifier 5 to drive each of the nine generating coils 3 sequentially The amplifier 

3S 5 outputs a 10kHz drive signal of 3amps mnis which causes the particular generating coil being driven to generate a 
quasi-static magnetic field. The frequency of the drive signal is chosen so that, within the range over which the locatbn 
of the sensor is to be detemnined, the field generated is a near-field electro-magnetic field i.e the wavelength is long 
compared to the distance from the generating coil to the sensing coll. 

Furthermore the drive signal frequency must be chosen so as to provide a compromise between sensor coil sen- 

^ sitivity, and the detrimental effects of electro-magnetic noise due to induced eddy currents within electrically conductive 
objects within the positioning range, since both of these aspects increase with frequency In the absence of electrically 
conducting objects a frequency of several hundred kilohertz may be used giving good sensor sensitivity and thus good 
range and positioning accuracy. In the presence of highly conductive objects, this frequency may need to be reduced 
to a few hertz. In this case a sensor coil may no longer be appropriate and may be replaced by an alternative magnetic 

^ field sensor, such as a flux gate magnetometer. In this embodiment a drive frequency of 10 kHz has been found to be 
a suitable compromise between sensitivity and immunity to interference from electrically conductive objects. 

Once the quasi-static field from a particular generating coil 3 is established, the value of the voltage Induced in 
the sensing coil 7 by this field is measured by the measurement unit 8. The signal from the sensing coil 7 is first amplified 
and then sampled at 40 kHz by a 16 bit analcgue-to-digital converter. The sampled signal is windowed using a Black- 
so man-Harris window, the 10 kHz component is extracted by the nnatched filter and hence a value representing the 
voltage induced in the sensing coil 7 is established. This value is passed to the controller 6 which stores the value and 
then instructs the amplifier 5 to stop driving the present generating coil 3 and to start driving the next generating coil 
3. When all nine generating coils 3 have been driven, or energised, and the corresponding nine voltages Induced in 
the sensing coil 7 have been measured and stored, the controller 6 calculates the locatbn and orientation of the sensor 

ss 7 relative to the field generators 1 and displays this on a display device 9. This calculation can be carried out while the 
subsequent set of nine measurements are being taken. Thus, by sequentially driving each of nine generating coils 3. 
arranged in three groups of three mutually orthogonal coils, the bcation and orientation of a single sensing coil 7 can 
be detemnined. 
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In order to describe the algorithm employed by the controller 6 to calculate the Icx^atbn and orientation of the 
sensor 7. a coordinate system will first be defined. In Figure 2 is shown a sensor, located at position P, whose axis is 
orientated along direction In general in order to determine the location and orientation of a single sensing coll within 
a field the x. y, z cartesian coordinates of the sensor and the elevation angle 9, and rotational angle 4>, must be found 
5 (see Figure 2). The vector distance R of the sensor from the origin. O, of the coordinate system is also shown in Figure 
2. Both the locatbn and orientation of the sensing coil within the field will affect the voltage induced in the coil by the 
field, but rotation of the coll about its axis will not affect the induced' voltage and thus does not constitute a further 
unknown quantity. 

Assuming now that a single field generating coil 3 Is placed at the origin O of the coordinate system with Its axis 
10 directed along the z-axis. When the generating coil is energised a field will be produced at the sensor location P which 
has a magnetic flux density B. With reference to Figure 3 this magnetic flux B can be resolved along the three axes of 
the coordinate system to give Bx, By and Bz and subsequently resolved along the axis of the sensor thus: 

B^ = B^cos(^ + BySin♦ (1) 



and 



so 



Bj = B J cos 6 + Bjjy Sin 6 



(2) 



25 



The voltage Vg, Induced in the sensor is related to the flux density via = Bg where is known and is a function 
of the frequency of the field and the characteristics of the sensing coil. It therefore follows from (1) and (2) that the 
voltage induced in the sensor at any x-y-z location and for any B-^ orientation is given by, 



Vg = kgtBj cos e + sin e(B^ cos + By sin $)) 



(3) 



30 



Formulae defining B^, By and B^ are developed from standard near field electronnagnetic theory In Appendix-A. 
Upon substituting the terms for B^,, By. B^ from equations (A-12) to (A-14) Into (3), it can be shown that. 



3S 



f2g2- - eosfl * 3g slnBfx eoBib ♦ v «ir>Al]-. (4) 



2vS/2 



3- 



where k^ is known and is a function of the cunent through, diameter of. and number of turns on the generating coil. 
The five unknown quantities, x, y, z, 8 and ^ are evident in (4): all other variables are known. 

Equation (4) has been derived, as stated above, for the case of a single generating coil 3 directed along the z- 
axis, there will of course be a corresponding equation for each of the three generating coils 3 of each of the three field 
generators 1 . 

It has been found that despite the complexity of this term (4) it is possible to determine the location and orientation 
of a single sensing coil by sequentially energising each generating coil. To explain this approach to position location 
the two dimensional case will first be considered. 

Figure 4 shows a field generator comprising two orthogonal field generating colls and Dy located at the origin 
of the coordinate system. The single sensing coil sensor is located at P and Its axis is parallel to direction S. The angle 
a is the angle between vector direction R of the sensor from the origin, and the direction S of the sensor axis. 

The voltages induced in the sensor when coils and Dy are energised sequentially are respectively, 

VfiDx = M^RDx • BeoxSina) (5) 

and 
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TO 



(6) 



where the and Dy sub-suffices relate to the field generated by the and Dy coils. Upon substituting .(A-1 ) and (A- 
2) from Appendix-A, (5) and (6) beconne, 

k k 

V^Thr = (2cose cosa - sine sina) (7) 
FT 

and 



k k * 
V«n« = (2sine cosa + cos8 sina) (8) 
IS Pi\ 

It has been noticed that the value of V(VsOx^ + V^Dy^) remains constant for a constant value of a. 
From (7) and (8) we can write. 

"/(VsDx^ + VgDy^) = ^ (4cos^ecos^a + sin^esin^a 
R 

-4sinesinacosBcosa . 
44sin^ecos^a +cos^Gsin^a 
446inGcosacosesina)** (9) 



which reduces to, 

j;;->-^=i^»777wa (10) 

35 

This can be thought of as corresponding physically to a circle of constant induced voltage in the sensor, centred 
on the field generator at the origin and lying in the x-y plane. This concept is shown schematically in Figure 5. If the 
t^pJndiidduaLmeasurements of induced voltage and measured at the spnsnr arfi us ed to calculate 
+ Vgjj^^ a circular or rotating field of constant strength can be simulated since ^V^^^ + represents the 
maximum voftage that could be induced in the sensor if a rotating field were used. This is desirable since equation 10 
gives a very simple relationship between R and a. 

The extension of this analysis to three dimensions is readily performed mathematically and conceptually very 
powerful because the approach taken does not require the axis of the generated field to be steered towards the sensor, 
but simply requires sequential energising of the individual generating coils. Thus for position determination in three 
^ dimensions of a single coil sensor, assuming three mutually perpendicular generating coils located at the origin of the 
coordinate system, we have 

so n 

It should be noted that the term TlTscos^a can only take values between 1 and 2. ignoring negative soluttons 
and thus any value of R computed from (11) is only weakly dependent on the value of a. For example, if a is assumed 
to be 7U2 whereas its correct value is zero, the value of-B-compuled from (11) is 80% of its c orrert vafuR . This in fact 
^ represents the worst case scenario since ou=0 means Vl + 3cos a = 2, while a=it/2 means 1 + 3cos a = 1 . 

Hence for each of three field generators a bounded value for R, the vector distance of the sensor from that particular 
field generator, can be calculated without any knowledge of the orientation a of the sensor. Since there are three field 
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generators located at different known positfons in the same plane (the x-y plane) and the distance R from each of them 
to the sensor has been calculated, the x-y-z coordinates of the sensor can be determined from simple trigonometry. 
This posittoning methodology is shown schematically in Figure 6. The three simulated spheres of constant induced 
voltage centred on each of the three field generators, and bounded by the potential error in R. overlap at two regbns. 

5 One regbn Is above the plane of the field generators and the other is below. In most applications, one solution Is clearly 
erroneous and the location of the sensor can easily be uniquely determined. 

At this stage the location of the sensor (but not its orientation) has been calculated to a limited degree of accuracy 
For some applications this may be adequate, but in general the precise location and probably orientation of the sensor 
are required. This is achieved by the use of an iterative procedure in which the estimate of the x-y-z coordinates of the 

10 sensor, derived from the values of R for each of the three field generators, are used In the appropriate equation (4) for 
each of the nine generating coils to estimate values of 9 and ^ for the sensor, from these a Is calculated for each of 
the three generators. Although 9 and ^ could be calculated from only two versions of equation (4), all nine versions 
are employed to improve the rate of convergence of the solution and its immunity from noise. The three values of a 
can then be invoked in the appropriate equation (11) for each field generator to calculate an improved estimate for R 

16 for each of the generators. This process Is repeated, progressively reducing the error in R and a for each generator 
until the desired levei of accuracy is achieved. It should be noted that this technique avoids the problems of non- 
convergence which would arise if equation (4) were utilised directly because a good estimate for R has been found 
before equation (4) is employed, and the estimate for R is bounded as shown schematically in Figure 6. 
In sumnnary, and with reference to Figure 7, the algorithm utilised by the controller 6 is as follows:- 

20 

1 . Assume a:=0 initially This ensures an over-estimate of R which guarantees an intersection of the radial distances 
from the three generator. 

p K/ioacnrQ iho x/r>itgrjoc induced in the sensor by each of the 9 indrvkJual generator coils, and then compute 
VVgQx^ + ^sDy^ + ^sDz^ ®ach Of the three generators. 
2S 3. Invoke a in (11 ) and compute R for each of the three generators. 

4. Compute the x-y-z coordinates of the sensor from the three values of R. 

5. Invoke these coordinates in the appropriate version of equation (4) for each of the nine generating coils and 
compute an improved estimate of 6 and This can be achieved by the use of. for example, the Gauss-Newton 
Least Squares optimisation technique. 

30 6. Use the improved estimates of 9 and ^ to cateulate a for each generator. 

7. Return to step 3 until the difference between the new and previous estimates of a reaches a sufficiently low 
value commensurate with the required positional accuracy in the x-y-z coordinates being achieved. 

Figure 8 depicts schematically the three spheres of constant induced voltage when the errors in R have been 
3S . reduced to altow the location of the sensor to be determined uniquely The technique employed thus guarantees con- 
vergence to a unk)ue location, with a precision that can be chosen in accordance with the requirements of the appli- 
cation. Indeed, it should be noted that In applications where the sensor is oKsving within the magnetic field, the number 
of iterations can be chosen dynamically for each calculation of the location of the sensor, thereby improving the effi- 
ciency of the process. For example, the first placement of the sensor typically requires 10 iterations before the solution 
40 is considered to have converged: this is considered to be so when the mean-square difference between the present 
and previous values of a is less than 10*®. Even with rapid movements of the sensor. It is unlikely that its angle a will 
change markedly from one positional placement to the next. By using the final value of a arrived at during the first 
placement as the initial estimate in the second placement, the number of iterations required to achieve the same 
convergence is significantly reduced. And so on for all subsequent placements. Experiments have shown that as few 
45 as 3-5 iterations are required for convergence after the initial placement. 

Although the algorithm described above with reference to Figure 7 ensures convergence to a unique location, 
allows both the location and orientation of a single coil sensor to be determined, and has proved to be robust even in 
the presence of noisy signals from the sensor coil 7, a second, alternative algorithm has been developed which has 
further advantages. 

so The first algorithm requires, at step 5, the solution of nine simultaneous equations relating 9 and ([> for each of the 

field generators to the estimate of the x. y and z coordinates of the sensor. This calculation can. dependent on the 
processing power of the controller 6, be time consuming, hence a second algorithm which is less computationally 
intensive has been developed. This algorithm enables the location and orientation of the sensor 7 to be determined 
more rapidly. The second algorithm Is based on the realisation that mathematically the voltages induced in the sensor 

ss 7 by each set of three generating coils 3 comprising each generator can be treated as vector quantities. This mathe- 
matical treatment enables an angle \j/ between the magnetic field lines and the directton vector of the sensor from a 
generator to be calculated. Once the values of v ^or each generator have been found there Is no need to employ 
equation (4) since the values of a can be catoulated directly from the values of y given a knowledge of the form of the 
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magnetic field. Since nine versions of equation (4) need no longer be solved this algorithm is coinputationally less 
intensive than the algorithm of Figure 7. 

The second algorithm will now be described in greater detail. In order to explain the algorithm clearly and to dem- 
onstrate the nnathematical insight on which it is based, the roles of the generating coils 3 and sensor coil 7 will be 

5 reversed i.e. for the purpose of the calculation the single axis field sensor 7 will be replaced by a continuously energised 
single axis field generating coil and the three orthogonal three-axis field generators will be replaced by three orthogonal 
three-axis field sensors. This is shown in Figure 9. Although it should be stressed that the reversal of roles here is 
simply for the purpose of mathennatical elegance, this reversed configuration will in practice be feasible and in some 
position location applications may be desirable. 

10 Referring now to Figure 9, let the vectors joining each three-axis sensor (10) to the single axis generator (11) be 
32 Ss angles between these vectors and the generator be c^, 02 and Og. The field produced by 

the single axis generator (11) will pass through each three-axis sensor (10), and the magnitude and direction of the 
field may be determined by processing the signals produced by each of the three orthogonal sensor coils (12). fonming 
each three-axis sensor (1 0). in response to the field. Let the signals in each of the threenaxis sensor (1 0) be represented 

16 by the vector quantities , and V^, where each component of the vectors corresponds to the signal in each of the 
orthogonal sensing coils (12). Let the angle between the field at each three-axis sensor (10) and the vectors £2 
and R3 be >|/^, and ^3 respectively, as shown in Figure 10. 

For the first estimate of the position of the generator 11 , the orientation of the generator (11 ) is unknown, and a, , 
ot2 and 03 are assumed to be zero. The magnitude of the vectors R^, and R2 and R3 are then calculated from equation 

20 (11). As for the first algorithm, because of the nature of equation (11) a bounded value for the distance of the generator 
(11) from each of the three-axis sensors (10) is found and the overlap of these bounded values can be used to give 
an initial estimate of the x, y and z components of each of the vectors R^, R2 and R3, 
The angles y-i, ^2 and ^3 are then calculated using the dot product, as follows: 



25 



30 



40 



45 



SB 



m 



Having found v„, we need to find ci„ to improve the estimate of position. Referring to Figure 11 y is the known 
angle and a is the required angle, d represents the calculated distance from the generator to the sensor. 

Since the generator is a simple dipole, the field at the sensor is given from equations (A-1 ) and (A-2) of the Appendix 

ss by: 



3, -|j^jco8a 



The angle of the field at the sensor is given by: 



so tanw = = 7i tana 

and so a is obtained from v using: 

tana„ =-2tanVn 

Having found a new estimate for a„. a new estimate of the generator posttbn is calculated using equation (11). 
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The process is repeated until the position converges to the required degree of accur&cy. 

Once the position of the generator (1 1 ) has been determined in terms of R„ and the orientation of the generator 
may be calculated in terms of 6 and <^ as follows. 

Let jJ be a unit vector defining the orientatbn of the generator relative to the sensors. Using the dot product, we 
can set up three equations to determine the three unknowns in U. 





R,. 


y = |2i| lyicosa, =|Rj| 


cos a. 


10 




U = |R^ iyicosa2 = |R2j 


COSOg 




R3. 


y = |R3j iyicosa3 = |Rgj 


COSttg 



IS 

These linear equations are solved to f ind U, and then the orientation In terms of 6 and ^ is given by: 



20 



if ^ U. 



25 



30 



^ » axctaxJ 0, ) 



35 



40 



(note that a four quadrant arctan functbn should be used). 

Although the formulation of the second algorithm has thus far been for the case of a single axis generator and 
multiple axis sensors the algorithm can be applied to the case of a single axis sensor and multiple axis generators. 
The only modification required between the two cases is the method by which the raw data for the algorithm (i.e. the 
voltages induced) is acquired. The equations developed above are directly applicable to the single axis sensor multiple 
axis generator case since the magnetic coupling between two coils is the same irrespective of which of the two coils 
is being driven. 

The steps to be followed when employing the algorithm for the single axis sensor and multiple axis generator case 
will now be summarised: - 



45 



1. Sequentially energise each of the three generator coils in each of the three generators 1 and measure the 
voltage induced in the sensor coil 7 by each generator coil i.e. measure V^^* V^y, V^^* ^Zx* ^2z* ^ax' ^Sz- 

2. Invoke On in equation (11) and compute |Rj for each of the generator 1. 2 and 3. (for initial estimate set a=0). 

3. From the Intersection of three spheres of radius |R J cateulate the vector quantities R^ , R2 and ^3. 

4. Taking the three voltages induced in the sensor coil 7 by a single generator 1 as a vector quantity e. g. ^ 



so 



calculate the angle of the field from the dot product V^. R^. 

5. Calculate the angles 0^ between the vectors R^ and the sensor axis from and equations A-1 and A-2. 
55 6. Repeat steps 2 to 5 until the desired level of positioning accuracy has been achieved. 

7. Use final values of and R^ to calculate the orientation of the sensor coil in tem^s of 6 and 

It has been found that use of the second algorithm can improve the speed with which the location and orientation 
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of a sensor is determined by a factor of approxinnately 15 compared to the first algorithm. 

For both algorithms the location and orientation of more than one sensor can be determined without the need to 
replicate the field generators 1 and amplifier 5. The field generated by any one field generating coil is measured at 
each of the sensors and the location and orientation of the sensors are simultaneous and independently calculated. 
5 The positions of the sensors may of course all be displayed on a single display unit 9. 

The simple, small sensor used in this embodiment means that it can provide position location in many situatbns 
where there is insufficient space for the three coil orthogonal sensor used in prbr art position location systems. A 
particular field of application is the medical field, where access through body vessels is required, for example in en- 
doscopy or non-invasive cardiovascular heart surgery. In these medical situatbns the present location system may 
10 replace the use of x-ray imaging (fluoroscopy), giving considerable advantages in cost and eliminating x-ray exposure 
to both patients and medical staff. The low frequency nragnetic fields used by the present system render the'human 
body transparent, while the use of low field strengths ensues the system Is intrinsically safe. 

During endoscopy it is desirable to know the path of the endoscope through the body This nnay be achieved using 
the present location system in three ways. Firstly the single sensing coil nnay be pulled along the biopsy tube and Its 
IS position at regular intervals along the tube stored and displayed to provide a 3D map of the path. Secondly, a tube 
containing approximately a dozen single coll sensors may be placed in the biopsy tube of the endoscope and the 
location of each of the sensors determined. This would be a retrofit to existing endoscopes. Allematlvely, the single 
coil sensors may be placed in the wall of the endoscope during manufacture. In the second two cases a real time 
picture of the path of the endoscope would be available at all times to the endoscopist. 
20 The present positioning system has been utilised in clinic field trials to Image in three dimensions the total config- 
uration of a colonoscope within the human abdomen. A sensor according b the present inventbn was placed inside 
the biopsy channel of an endoscope. 

The small inner diameter of the biopsy channel, typically 3. 7mm for a colonoscope, not only dictates that the 
sensor be of vanishingly small diameter, but also that it may only comprise a single coil, typically 1cm in length, orien- 
ts tated along the axis of the instrument. The algorithrns of the present positioning system processes the signals from 
this sensor in such a way as to calculate the position of the sensor within the biopsy channel independent of its orien- 
tation. Such independence is cmcial in colonoscopy since the sensor may adopt any orientation for a single x-y-z 
location. 

The positioning algorithm resides as software within an IBM 486 personal computer which, upon processing the 

30 information taken from the sensor at numerous discrete positions along the biopsy channel, then displays the path 
followed by the sensor as a continuous line on the nnonitor. Clearly this path corresponds precisely to that of the en- 
doscope. Moreover, because the information from the sensor at each location relates to three dimensbns, the imaged 
path on the monitor is likewise displayed in three dimensions. Visually the system achieves this by the use of 'grey 
scale" colour coding whereby portions of the path further from the viewer (i. e. down into the screen) appear in darker 

3S shades of grey than the "under* portion. This feature is unique among all conventional imaging techniques for colon- 
oscopy and represents a major advance in the field. 

To display the path of the endoscope, the endoscopist first passes the sensor ctown the biopsy channel until it 
reaches the tip of the endoscope. For convenience we have encapsulated the sensor within a hollow tubular catheter 
of the type used routinely with endoscopes. The catheter is then withdrawn at a uniform speed (although this is not 

40 critical) while the system repeatedly determines the position of the sensor at discrete instances during its motion. During 
withdrawal the path of the Instrument is displayed on the monitor in three dimensions. In many situations a total image 
of the endoscope is not required, In which case the sensor need only be withdrawn along that portion of the instrument 
of interest. To cater for patients lying in a variety of positions, perhaps changing during the investigation, the image 
may be rotated in any direction. This is particularly advantageous in establishing the radius of curvature of any bend 

45 in the endoscope that happens to lie along the viewing axis. For example, a bend that is In fact gradual, and hence 
poses no concern, can appear abrupt if viewed from some directions. A useful zoom facility on the image is also 
provided. When the system \s in normal use, the system display would ideally be sited next to a standard camera 
monitor used to display the view from the endoscope. In this way the endoscopist Is conveniently presented with the 
path of the instrument in three dimensions on one display, and the internal view from the endoscope optics on the other. 

so Initial validation of the system was performed with the aid of a rigid plastic framework to hold the endoscope in 

one of a number of predefined configuratbns. X-ray imaging and the present magnetic field system were applied to 
seven different configurations of the endoscope. These Included a sigmoid loop, an alpha loop, a reverse alpha loop, 
a gamma bop. and an "N" loop. The results, three of which can be seen in Figure 12 showed close agreement between 
the image produced by the present positioning system (shown on the left) and the X-ray innage (shown on the right) in 

ss each case. The nature of the overlapping portbns of the colonoscope can be clearly seen from the images produced 
by the present positioning system. Some distortion of the images was caused by the metallic construction of the colon- 
oscope perturbing the magnetb fields. However, this was minimal and the cobnoscope configuration Is clearly evident 
from the images.' 
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The clinical trails involved three patients undergoing colonoscopy for a number of different indicatbns. Ethical 
approval was obtained, as was written consent. The patients were sedated with a combination of pethidine and mida- 
zolam before the examination. The colonoscope used was a Pentax type FC38LH. 

For the majority of each examination, the sensor was fully inserted into the biopsy channel, and the display was 

s configured to show the progress of the tip of the endoscope in real time. When progress became difficult, the sensor 
was withdrawn, which immediately produced an image on the screen of the total path of the endoscope. With the aid 
of this image the removal of bops was straightforward, by using clockwise or anti-clockwise twist and simultaneous 
withdrawal of the endoscope. Similarly, when re-inserting the instrument the refonnation of loops was prevented by a 
combination of abdominal pressure and torque. Where abdominal pressure was required the sensor was positioned 

10 inthe loop, so enabling the endoscopist to see, by referring to the displayed image, whether pressure was being applied 
in the desired direction and to the correct extent. In each case examination around to the caecum was achieved (i.e. 
total colonoscopy) and the procedure was tolerated well by the patients. During the examinatbns, X-ray pictures were 
taken for comparison against those obtained with the nr^gnetic system. Two of these, a plan and side view, are shown 
in Figure 13 together with the corresponding image from the magnate system. Agreement between the two is very 

IS closB, the devlatk>n being largely attributable to patient movement between the two exposures. 

The system has been shown to image the configuration of the endoscope within the patients's abdomen with close 
agreement to the X-ray image. The three dimensionality of the Image has proven to be of great help in deciding the 
strategy for removing loops which form in the path of the endoscope during intubation. Indeed, this improvement in 
visualisation is likely to be of great benefit in teaching colonoscopy, as well as enabling experienced endoscopists to 

20 improve their technique when facing difficult cases. The intrinsically safe nature of the system allows it to be in con- 
tinuous use throughout the examination, presenting the endoscopist with as many images as the circumstances require. 
. This contrasts markedly with fluoroscopy which can only offer images intermittently and carries an exposure time limit 
for reasons of patient safety, and X-ray pictures which are essentially only a "one-shot" option. Moreover, protective 
clothing need not be worn by any of those present at the examination while the system is in use, nor is it necessary 
for the examination room to be in any way specially prepared. Indeed, the system frees such examinations from having 
to take place in a room apart from the ward. If need be such examinations could be carried out in complete safety and 
with no loss In overall Integrity, at the patient's own bed in the ward. 

A number of medical studies have considered the efficacy of colonoscopy as a screening methodology in asymp- 
tomatic subjects and have shown a significant detection rate for adenomas and carcinoma in subjects over the age of 

30 60. Of particular note here is that some 50% of lesions were proximal to the splenic flexure, hence the importance of 
performing a total colonoscopy In such cases. The ability to conduct total colonoscopes routinely and efficiently is 
therefore an important objective. On the other hand it must be remembered that colonoscopy (total or othenvise) Is 
associated with a certain morbidity and mortality due to the need to apply mechanical stress during intubation or with- 
drawal. The overall Improvement in visualisatbn that the present system affords, particularly it's three dimensbnality, 

35 should both raise the efficacy of total colonoscopy and reduce the risk of perforation. This in turn may also help to 
reduce the dosage of analgesic and sedative drugs required. 

Although the apptk:atlon of the present positioning system to colonoscopy has been specifically addressed, the 
scope of the medical applications extends tar beyond this by virtue of the very small size of the sensor(s). For example, 
bronchoscopy, gastroscopy and procedures involving a nasogastric or endotracheal tube could all utilise the sensor 

40 described herein its present catheter form. Numerous other medical applications requiring position or orientatkan In- 
formation could benefit from either a single or multiple sensor implementation of the system. 

Data gloves which facilitate the location of a wearer's hands, are used in both medical and virtual reality applica- 
tions. They enable the position and direction of each of the fingers to be determined. The prior art magnetic field location 
system using a three coil orthogonal sensor Is clearly not applicable, so current data gloves use fibre optic strain 

45 gauges. These require calibration every 2-3 minutes. The ability to locate single coil sensors means that the sensors 
may be wound around each joint of each finger giving a system which is less bulky, more accurate and only requires 
calbratbn during the manufacture of the gloves. 

A particular area of application for the present positioning system comprises that of the so called "man-rnachine 
interface". There are numerous situatioris in which a hunr^an operator needs to interact with a machine, or computer, 

so normally comprising some form of display device, examples of such interactions are with a conventional personal 
computer, a video conferencing system, or a virtual reality environment in which the operators field of view is filled by 
the display device, which in this case may be three dimensional. The present positioning system allows an operator 
to wear small, single coil sensors about his body to enable his movements to be detected and interpreted by a machine 
without the need for physical contact between the operator and the machine. For example the positioning system of 

ss the present inventbn could enable an operator to interact with images on a televisbn or computer screen without the 
use of a conventional keyboard, mouse or stylus. The operator could wear single coil sensors on his fingertips, for 
example in thimbles, or a thin gkTve, the bcation and orientation of whbh could be detected within a magnetic field 
generated within the vicinity of the display screen. Linking the positioning system to the computing system would allow 
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the computing system to have knowledge of the position of the operators fingertips ih three dimensions. A computed 
drawn replica of the user's hand which precisely emulates the movements of the user's own fingers, could then be 
utilised by the user to interact with the computer system. Thus when the user makes hand movements the virtual hand 
on the screen can be made to grasp and manipulate objects in the display, for example moving portions of text, rotating 

5 an engineering drawing, selecting an icon to activate a software program, etc. The virtual hand could also be used to 
control windows and menus and to draw diagrams. The advantage of such a man machine interface is that its use is 
completely intuitive, requiring no training. 

Since the positioning system of the present invention enables the position of a sensor to be located In three di- 
mensions, the extension of such a man machine Interface to a three dimensional virtual reality environment is clearly 

10 possible. In this case the computer system involved may need information regarding the position of other parts of the 
operator* body than his hands, for example the image displayed to the operator may be dependent on the location and 
orientation of his head, in which case small single coil sensors can clearly be worn for example on each temple, 

In a second embodiment the invention enables a sensor, comprising three orthogonal sensing coils, to be locatei^ 
in three dimensions relative to a single field generator comprising three orthogonal generating coils. 

IS With reference to Figure 14, a fleki generator 1, comprising three generating colls 3, as prevk)usly described Is 
mounted on a surface 2. Each generating coll is electrically connected to an amplifier 5 and is driven as previously 
described. 

The sensor 7 in this embodiment comprises three mutually orthogonal sensing coils. A, B and C, each of which is 
separately electrically connected to a measurement unit 8. 

so In use the three generating coils are sequentially energised as previously described, but when each coil Is energised 
the voltages induced in each of the three sensing coils V^, Vg and are measured by the measurement unit 8 and 
stored by the controller 6. The controller 6 then calculates from these three voltages the location of the sensor 7 relative 
to the single field generator 1 . 

The controller is able to calculate the location of the sensor, even though the axes of the generated fields have 

2S not been directed towards the sensor, by employing an algorithm which weights the voltages induced in the three 
sensing coils by a kx:ation dependent weighting, and then alters these weightings to achieve a calculated maximum 
field strength at the sensor. In order to more fully describe this algorithm the field from a single small coil is first con- 
sidered. 

The nr^gnetic field produced by a small coil, from equations (A-1 ) and (A-2), is given by: 

30 

B = \ (2a^ cose + a^ sin 6) (12) 



3S where 

R rr distance from the coil 
9 = angle from the axis of the coil 
k = constant for coll (size, drive current, no. turns etc). 
40 Qr is a unit vector in the direction of Br (see Appendix A - Fig. A-1 ) 

Bq is a unit vector In the direction of Bq (see Appendix A - Fig. A-1) 

Now, the magnitude of the magnetic field 

^ |B|=^3V(3cos^e + 1) (13) 

R 

and so it can be seen that for a given distance from the coil, the field strength is greatest when 9=0 i.e. on the axis of 
so the coil. Clearly, ff the effective axis of the coil could be directed towards the sensor, the sensor would experience a 
maximum in field strength. 

In order to steer the effective axis of the coil without physically moving it. addittonal coils are required. To steer the 
effective axis over 3D, three coils are required In total. Assuming three mutually perpendicular coils D,j. Dy. lying 
along each of the cartesian axes x, y and z, each coil being centred on the origin, by setting the currents to each coil as: 

55 

\^ = I cos9 cos^ ^ 
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I = I cose sm^ 



I J = I sine 

s 

the effective axis of the resulting field may be steered without changing the nriagnitude of the field. ^ Is the angle 
anticlockwise from x in the xy plane, and 6 Is the elevation towards the z axis. 

Assuming the notation of Figure 2, OP represents the effective axis of the field. That is a single drive coil, centred 
10 on the origin, with Its axis along OP. fed with current I. would create the same field as the three coil arrangement with 
the currents 1^, ly and 1^ as described. 

Thus if the field strength at the point we wished to locate could be measured, we would find that when axis OP 
pointed at this point, the field strength would be a maximum. 

The field strength is measured using 3 orthogonal sense coils, centred on a single point. In this case an AC field 
IS must be used In the drive coils. Let the sensor coits be A, B and C, and let the amplitude of the voltages induced be 
V^, Vb and V^. The field strength can be computed from 

20 ' 

where kg = a constant for the sensor and frequency used. 

The effective axis of the resulting field, could be physically steered towards the sensor, and V;^. Vg, monitored 
to maximise B. However this is difficult in practice to achieve since both 9 and $ would need to be simultaneously 
altered while measurements from the sensor are taken. This leads to slow position location, and limits the system to 
2S locating a single sensor. The approach adopted In this embodiment is as follows. The drive currents for all the coils 
are set to I. and not to the values which would be required to physically steer the effective field axis, as discussed above. 
I.e. 

I, = l 

30 ly=| 
l,= l 

Effectively steering of the field axis is carried out AFTER the field measurements have been made by weighting 
or scaling these measurements by location dependent weighting factors. Thus, instead of physically altering 9, ^ and 
35 then measuring B, the following technique is used. 

1 . Switch on D^, with 1,^ = I 

2. [Measure Vadx. Vbdx. Vcdx 

3. Switch off Dxi Switch on Dy, with ly = I 
40 4. Measure V^Dy. Veoy. Vcoy 

5. Switch off Dyi Switch on D^, with 1^ = I 
a Measure V^q^, Vgo^, V^dz 
7. Switch off Dj 

45 For the physically steered field: 1^ = I cosecos(t), rather than I. The same result is achieved by weighting the results 

from step 3 by cosecos(t>. The same logic applies to the remaining results, using the relevant weighting factor. 
Thus: 

so = <<8^(((Vadx cos<fr + V^Dy sin<|))cose + V^q, sine)^ 

+ i(^BDK cos* + VgDy 6ln(|>)cose + Vbd^ sinsf 
+ ((VcDx cos* + VcDy sin*)cose + VcD^ sine)^) 

55 

Note that the "signs" of the amplitude are important 

ag. 
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f^ase shift = O +V6 



^ phase shift =n -ve 

In this expression for B^, e and $ are the only variables. 

in order to find the values of 6 and ^ which give the maximum B^, the Gauss-Newton optimisation technique is 
used. This copes well with sum of squares type expressions. The expression for is well behaved, and only a few 
10 iterations are required. 

In order to find the precise location of the sensor we must now find R. 

If. we square and sum the field magnitudes at the sensor for each generator coil, we find that: 

and so R may be found from: 

The cartesian coordinates of the sensor are then 

25 

X = R cosecos(|> 



y = R cos6sln<|) 



2 = R sine 

35 As with the first embodiment the location of multiple sensors Is facilitated because the generating coils are only 

energised sequentially allowing the generated field to be simultaneous measured at any number of Icications. 

Although in both embodiments of the invention described herein the voltages induced in the sensor coil 7 by the 
generating coils 3 are distinguished one from the other by employing a time multiplexing approach, i.e. the generating 
coils are energised sequentially, a frequency multiplexing approach nrtay also be adopted within the scope of the present 

40 invention. For example in such an approach each generator coil 3 could be driven at a different frequency so that a 
plurality of generating coils 3 could be simultaneously energised white still allowing the voltage induced in the sensor 
7 by each generating coil to be distinguished by its frequency. In such an arrangement the sensor would need to be 
responsive to all the energising frequencies and some form of frequency filtering would need to be provided. This 
filtering could be provided by discrete physical bandpass filters electrically connected to the sensor 7, or, if an A to D 

45 converter is employed as described herein, filtering of the signal from the sensor 7 can be accomplished by signal 
processing software in the controller 6. The use of frequency multiplexing to acquire the data for position determination 
can significantly increase the operating speed of the positioning system since measurements from generating coils 
can be taken simultaneously. Disadvantages of such a frequency multiplexing system are that it is more complex than 
a time multiplexed system and requires greater electrical bandwidth. A combination of time and frequency multiplexing 

so could of course be used. 

In both emt>odiments it is desirable that the quasi-static magnetic field generated by a coil is established quickly 
and is allowed to decay qutekly For this reason it is prefen^ed to use a first order rather than a second order drive 
circuit. For the generating coiis employed the field settles within one cycle of being switched on. 

it will be appreciated that alternative configurations of both embodiments for different applications, for example 

55 locating a sensor within a two dimensional plane, are envisaged within the scope of the present invention. 

As will be clear to one skilled in this art. the roles of the generating and sensing coils may be reversed while still 
benefitting from the advantages of the present Invention. That is the sensing coil or coils may be used as field generating 
elements, and the generating coils may be used as field sensing elements. 
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This reversal of roles has particular advantage where a static field, such as that generated by a bar magnet \t 
employed according to the first aspect of the invention, since such a field generating element must be effectively 
permanently "energised". The reversal of roles albws the "sensor" permanently to generate a field which Is sensed at 
each of the "generating elements" and the location and orientation of the "sensor" is then determined as before. 

Appendix A 

Consider a cunent. 1, flowing through a small planar coll of radius, b, (Fig A-1). The frequency of I is chosen to be 
sufficiently bw such that static field distributions apply. 




FIq A-1 . Resolving magnetic flux density at a distance. 

For a point P in the field whose distance R from the coil is such that R»b, it is readily shown for example in D K 
Cheng, Field and Wave Electromagnetics. 2nd Ed, Addison Wesley, 1989, that 



30 



Br = - 



2 cos e 



(A-1) 



3S 



and, 



sin 8 



(A-2) 



40 



4S 



where k^ is a known function of I and b, Br and represent the vector components of the magnetfe flux densityat 
point P resolved along axes parallel to the line R and the angle 6. Note that, by convention, 8 is measured from the 
axis of the coil. 

In order to resolve the magnetic flux density at P onto a 3-dimensional cartesian coordinate system, consider first 
the coil in the y-x plane, centred on the origin (Fig A-2). 



so 



ss 
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10 



IS 



20 




FiQ A-2. Resolving magnetic tiux density onto an x-v-z cartesian c oordinate system. 

It pobit P is distance xlronUhe coll (i.e. origin) along the x-axis and its vector dis tance is R. the distance in the y- 
z plane by Pythagoras is V?-?. Since = x2+ y2 + z2, this distance reduces to Vy + z .as shown. It then follows 
that. 



(A-3) 



ss and, 



30 Resolving the magnetic flux density at P onto a cartesian system gives, for the x-axis component, 



SB From (A-1 ) and {A-2) this becomes, 



40 



45 



and from (A-3) and (A-4), 



Bjj = Bp cos e - Bg sin 6 



B, = -4(2cos^ e-sin^ 8) 
R^ 



B, = ^5(2x'-y'.z^) 



Resolving similarly onto the y-z plane gives. 

B„ = Bo sin + B« cos (t> 

which from (A-1) and (A-2) becomes, 



(A-4) 



(A-5) 



ss 



B„ = -|(3 cose sin 6) 
^ R^ 



and from (A-3) and (A-4), 
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B„ = %(3x7y^l^) 



(A-6) 




and, 



30 



3S From (A-6) these become, 



B =^^(3xy) (A-7) 

40 • 

and, 

Bj^ = ^^(3xz) (A-8) 

45 ^ 

For a coil (dipole) in the y-z plane, equations (A-6). (A-7) and (A-8) f uffy describe the resolved cartesian components 
of the magnetic flux density at a point P located at a radial distance R from the coil. The corresponding equatfons for 
coils in the x-y and x-z planes can be devetoped in an identical manner. The complete set of fomiulae can therefore 
so be summarised thus. * 



For a coil in the y-z plane : 



55 



(A-9) 
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For a coil in the x-y plane > 

TO 



IS 



For a coil in the x-z plane > 

25 



By = ^(3xy) (A-IO) 
R 



B, = ^|3X2) (A.11) 
R 



B, = ^^(3X2) (A-12) 



By = -^(3y2) (A-13) 
R 

20 B, = ^s(22^-x^-y^) (A-U) 

' R^ 



9 



B, = -|(3xy) (A.15) 
R^ 



50 B„ = -!^(2y^-x^.z^) ^ (A-16) 

B, = ^(3yz) (A-17) 

35 R 



Claims 

40 1 . A method of determining the location and/or orientation of a field sensor (7) relative to a plurality of field generators 
(1) of known location, each flekJ generator (1) comprising a plurality of collocated field generating elements (3), 
each generating element (3) being energised and a measurement being made at the field sensor (7) of the re- 
spective field generated by the generating element (3). the method being characterised by the steps of: 

45 1 ) ) for each field generator (1) calculating, from the measurements of the field generated by each of its gen- 

erating elements (3) and an estimate of the orientatton of the sensor (7), an estimate of the distance from that 
particular field generator (1 ) to the sensor (7). and 

2) utilising the estimates of distance from step 1 ), and the known location of the tieki generators (1 ) to calculate 
an estimate of the location of the sensor (7); 
50 3) employing the estimated location of the sensor (7) from step 2) and the measurements of the field at the 

sensor (7) to calculate a new estimate of the orientatbn of the sensor (7); 

4) repeating steps 1 ) to 3) iteratively, with step 1 ) enriploying the new estimate of sensor orientation from the 
preceding step 3), to improve the estimates of location and orientation of the sensor (7). 

55 2. A method as claimed in claim 1 . wherein each generating element (3) of each field generator (1) is energfeed in 
sequence so that only one generating element (3) is energised at any given time. 

3. A method as claimed In claim 1 . wherein at least two generating elements (3) are simultaneously energised, each 
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of these at least two generating elements (3) being energised at a different frequency from the other simultaneously 
energised generating element(s) (3). 

4. A method as claimed in claim ^ wherein at step 3). the estimated location of the sensor (7) calculated in step 2) 
5 and the measurements of the field at the sensor (7) are used to calculate, for each field generator (1 ), the direction 

of the respective field at the sensor (7) and from these field directions calculating a new estimate of the orientation 
of the sensor (7). 

5. A method as claimed in claim 1 or claim 4 wherein, at step 1 ), the initial estimate of the orientation of the sensor 
10 (7) is taken to be the orientation as calculated by the immediately preceding application of the method or, If this is 

the first application of the method, the axis of the sensor (7) is initially assumed to be directed towards the respective 
field generator (1). 

6. A method of determining the location of a field sensor (7), comprising a plurality of collocated field sensing elements 
IS (A,B.C). relative to a field generator (1 ) comprising a plurafity of collocated field generating elements (3), each field 

generating element (3) being energised in turn and the value of the resultant field strength being measured at the 
field sensor (7), the measured field strength being dependent on the location of the sensor (7) within the respective 
field, the method being characterised by the steps of: 

20 1) calculating, by utilising all the values of measured field strength and an estimate of the direction of the 

sensor (7) from the field generator (1 ), a direction dependent weighting factor for each field generating element 
(3) so that the calculated field strength B is equal to the field strength B that would exist at the sensor (7) if 
the axis of the field were directed towards the sensor (7), 

2) iteratively altering the direction dependent weighting factors to maximise B and thus to determine to a 
2S desired level of accuracy the direction of the sensor (7) from the field generator (1 ). and 

3) employing the measured values of the field strength to calculate the distance of the sensor (7) from the field 
generator (1) and hence, from the direction of the sensor (7) in step 2). the location of the sensor (7) relative 
to the field generator (1 ). 

30 7. A method of determining the location and/or orientation of a field generator (7) relative to a plurality of field sensors 
(1) of known location, each fieki sensor (1) comprising a plurality of colkxated field sensor elements (3), the gen- 
erator (7) being energised and a measurement being made at the field sensor elements (3) of the respective field 
generated by the generator (7), the method being characterised by the steps of: 

35 1) for each field sensor (1) calculating, from the measurements of the field generated at each of its sensor 

elements (3) and an estimate of the orientation of the generator (7), an estimate of the distance from that 
particular sensor (1) to the generator (7), and 

2) utilising the estimates of distance from step 1 ), and the known locatfon of the field sensors (1 ) to calculate 
an estimate of the location of the generator (7); 
40 3) employing the estimated location of the generator (7) from step 2) and the measurements of the field at the 

sensors (1) to calculate a new estimate of the orientation of the generator (7); 

4) repeating steps 1) to 3) iteratively, with step 1) employing the new estimate of generator orientation from 
the preceding step 3), to improve the estimates of locatbn and orientation of the generator (7). 

^ 8. Apparatus for determining the location and/or orientation of a field sensor (7) comprising: 

a plurality of field generators (1 ) each comprising a plurality of collocated field generating elements (3), 
energising means (5) operable tor each field generator (1) to energise each of the plurality of field generating 
elements (3) to produce a field within which the location and orientation of the sensor (7) may be detemnined, 
so measurement means (8), connected to the field sensor (7), to measure and output a parameter of the field at 

the sensor (7), and 

control means (6) to control the energising means (5), the control means (6) connprising storing means to store 
parameters output by the measurement means (8), and calculating means to calculate, from the stored meas- 
urements of the fields, the locatbn and orientation of the sensor (7), 
ss the apparatus being characterised in that the calculating means comprises: 

bcatlon estimating means being operable to calculate, from the stored measurements of the fields and 
an estimate of the orientation of the sensor (7), an estimate of the distance of the sensor (7) from each 
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o1 the field generators (1) and. using the estimates of distance and the khown location of the field gener- 
ators (1), to calculate an estimate of the location of the sensor (7); 

orientation estiniating means being operable to calculate, from the estimate of sensor location by the 
location estimating means and the stored measurements of the field, a new estimate of the orientation of 
5 the sensor (7); 

means to iteratively improve the estimate of the location and orientation of the sensor (7), which means 
in use repeatedly employ the location and orientation estimating means in sequence, wherein the output 
of the orientation estimating means is adopted by the location estimating means as the new estimate of 
sensor orientation at each iteration. 

10 

9. Apparatus for determining the locatbn of a field sensor (7) as claimed in claim 8. wherein the field sensor (7) 
comprises a single field sensing element and the parameter of the field sensed by said single field sensing element 
Is dependent on the orientation of the field sensing element within the field 

T5 10. Apparatus for determining the location of a field sensor (7) as claimed in claim 8 or claim 9, having three field 
generators (1 ). each connprising three mutually orthogonal collocated field generating elements (3). 

1 1 . Apparatus for detemnining the location of an object relative to a field generator (1 ). the apparatus comprising 

20 a field generator (1 ) comprising a plurality of field generating elements (3), 

energising means (5) to sequentially energise each of the field generating elements (3) to produce a field 
within which the object may be located, 

a sensor (7), to be attached to the object, comprising a plurality of collocated, field sensing elements (A,B,C). 
each able to sense a parameter of the field, i 
25 measurement means (8) to measure and output the parameter of the field sensed by each field sensing element 

{A,B,C) when each field generating element (3) is energised, and 

control means (6) to control the energising means (5). store the parameters output by the measurement means 
(8) and calculate the location of the sensor (7). relative to the field generator (1). 
the apparatus being characterised in that the control means (6) comprises 

30 weighting means to weight each of the parameters, sensed by the sensor (7), by a direction dependent weight- 

ing factor for each field generating element (3). so that the calculated field strength B is equal to the field 
strength B that would exist at the sensor (7) if the axis of the field were directed towards the sensor (7). 
means to iteratively alter the direction dependent weighting factors until a maximum value of the field strength 
is found, and to determine the direction of the sensor (7) from the field generator (1) to a desired level of 

3$ accuracy, and 

calculating means to calculate the distance of the sensor (7) from the generator (1 ) and hence the location of 
the sensor (7) relative to the generator (1 ). 

1 2. Apparatus for determining the location of an object as claimed in claim 11 , wherein the field generator (1 ) comprises 
40 three mutually orthogonal field generating elements (3) and the sensor (7) comprises three mutually orthogonal 

field sensing elements (A.B,C). 

13. Apparatus as claimed In any one of claims 8 to 12, arranged for use in surgery, therapy or diagnosis of a human 
or animal body. 

45 

14. An endoscopy system for use in the endoscopy of a human or animal body, the endoscopy system comprising 
apparatus as claimed In any one of claims 8 to 12. 

15. Apparatus for monitoring the position of an endoscope within a human or animal body, the apparatus comprising 
so a field sensor (7), 

a plurality of field generators (1) each comprising a plurality of field generating elements (3), 
energising means (5) operable for each field generator (1) to energise each of the plurality of field generating 
elements (3) to produce a field within which the location of the endoscope may be determined, 
55 ' measurement means (8), connected to the field sensor (7). to measure and output a parameter of the field at 

the sensor (7). and 

control means (6) to control the energising means (5), store parameters output by the measurement means 
(8), and calculate the location of the sensor (7) relative to the plurality of field generators (1), the apparatus 
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being characterised in that the control means (6) are arranged to calculate the'locatbn of the sensor (7) relative 
to the plurality of field generators (1 ) by the method as claimed in any one of claims 1 to 7. 

16. Apparatus lor determining the location and/or orientation of a sensor (7) as claimed in any one of claims 8, 9 or 
5 1 0, wherein one or more of said sensors (7) are attached to parts of a glove, the output of said apparatus providing 

data defining the location and orientation of the respective parts of a glove wearer's hand. 

17. Apparatus as claimed in claim 16. wherein one or more sensors (7) are attached to fingers of the glove. 

10 

Patentanspruche 

1. Verfahren zum Bestimmen des Ortes und/oder der Orientierung eines Feldsensors (7) relativ zu mehreren Fel0- 
erzeugern (1) an bekannten Orten, wobei jeder Felderzeuger (1) mehrere zusammengestellte feiderzeugende 
16 Elements (3) enthatt. wobei jedes erzeugende Element (3) eingeschaltet wird und bei dem Feldsensor (7) eine 
Messung des jeweiligen durch das erzeugende Element (3) erzeugten Feldes ausgef uhrt wird, wobei das Nferfahren 
gekennzeichnet ist durch die folgenden Schritte: 

1) fur jeden Felderzeuger (1) Berechnen einer Abschatzung der Entfemung von jenem besonderen Felder- 
20 zeuger (1 ) zu dem Sensor (7) aus den Messungen des durch jedes seiner erzeugenden Elemente (3) erzeug- 
ten Feldes und einer Abschatzung der Orlentiemng des Sensors (7). und 

2) VenAfenden der Entfernungsabschatzungen aus Schritt 1) und des bekannten Ortes des Felderzeugers (1) 
zum Berechnen einer Abschatzung des Ortes des Sensors (7); 

3) Ven^/enden des abgeschatzten Ortes des Sensors (7) aus Schritt 2) und der Messungen des Feldes bei 
2S dem Sensor (7) zum Berechnen einer neuen Abschatzung der Orientierung des Sensors (7); 

4) iteratives Wiederholen der Schritte 1) bis 3) unter Venwendung der neuen Abschatzung der Sensoforien- 
tierung aus dem vorausgehenden Schritt 3) In Schritt 1) zum Verbessem der Abschgtzungen des Ortes und 
der Orientierung des Sensors (7). 

30 2. Verfahren nach Anspruch 1, bei dem jedes erzeugende Element (3) jedes Felderzeugers (1) sequentiell einge- 
schaltet wird, so da8 zu irgendeiner gegebenen Zeit nur ein erzeugendes Element (3) eingeschaltet ist. 

3. Verfahren nach Anspruch 1 , bei dem zumindest zwei erzeugende Elemente (3) gleichzeitig eingeschaltet werden, 
wobei jedes dieser zwei erzeugenden Elemente (3) bei einer von dem (den) anderen gleichzeitig eingeschalteten 

3S erzeugenden Element(en) (3) verschiedenen Frequenz eingeschaltet wird. 

4. Verfahren nach Anspruch 1 . bet dem im Schritt 3) der In Schritt 2) berechnete geschatzte Ort des Sensors (7)'Und 
die Messungen des Feldes bei dem Sensor (7) venwendet werden, um fur jeden Felderzeuger (1) die RIchtung 
des jeweiligen Feldes bei dem Sensor (7) zu berechnen und aus diesen Feldrtchtungen eine neue Abschatzung 

40 der Orlentiemng des Sensors (7) zu berechnen. 

5. Verfahren nach Anspmch 1 oder Anspruch 4, bei dem bei Schritt 1 ) als Anfangsabschatzung der Orientierung des 
Sensors (7) die durch die unmittelbar vorausgehende Anwendung des Verfahrens berechnete Orientierung ge- 
nommen wird oder, falls dies die erste Anwendung des Verfahrens ist. die Achse des Sensors (7) anfangs als zum 

45 entsprechenden Fekjerzeuger (1 ) gerichtet angenommen wird. 

6. Verfahren zum Bestimmen des Ortes eines Feldsensors (7) mit mehreren relativ zu einem Felderzeuger (1) zu- 
sammengestellten Feldsensorelementen (A. B, C) mit mehreren zusammengesteitten feklerzeugenden Elementen 
(3). wobei jedes feiderzeugende Element (3) der Reihe nach eingeschaltet wird und der Wert der resultierenden 

so Feldstarke bei dem Feldsensor (7) gemessen wird, wobei die gemessene Feldstarke von dem Ort des Sensors 

(7) innerhalb des jeweiligen Feldes abhangt. wobei das Verfahren gekennzeichnet ist durch die folgenden Schritte: 

1) Berechnen eines richtungsabhangigen Gewlchtungsfaktors fDr jedes feiderzeugende Element (3), soda8 
die berechnete Feldstarke B gieich der Feldstarke B ist, die bei dem Sensor (7) vorhanden ware, wenn die 

ss Achse des Feldes zum Sensor (7) gerichtet ware, unter VenAfendung aller Werte der gemessenen Fekistarke 

und einer Abschatzung der Richtung des Sensors (7) von dem Felderzeuger (1 ), 

2) Iteratives Andern der rk:htungsabhangigen Gewchtungsfaktoren. um B zu nnaximieren und daher die Rich- 
tung des Sensors (7) von dem Felderzeuger (1 ) bis zu einem gewunschten Genauigkeitsgrad zu bestimmen, 
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und 

3) Verwenden der gemessenen Werte der Feldstarke zum Berechnen der Entfemung des Sensors (7) von 
dem Felderzeuger (1) und daher aus der Richtung des Sensors (7) In Schritt 2) des Ortes des Sensors (7) 
relativ zu dem Felderzeuger (1). 

Verfahren zum Bestimmen des Ortes und/oder der Orientierung eines Fetderzeugers (7) relativ zu mehreren Feld- 
sensoren (1) an bekannten Orten, wobei jeder Feldsensor (1 ) mehrere zusammengestellte Feldsensorelemente 
(3) enthaft, wobei der Generator (7) eingeschaltet wird und wobel eine Messung des jeweiligen durch den Erzeuger 
(7) erzeugten Feldes bei den Feldsensorelementen (3) ausgefuhrt wird, wobei das Verfahren gekennzeichnet ist 
durch die tolgenden Schritte: 

1 ) f Qr jeden Feldsensor (1 ) Berechnen einer Abschatzung der Entfemung von jenem besonderen Sensor (1 ) 
zu dem Erzeuger (7) aus den Messungen des bei jedem seiner Sensorelemente (3) erzeugten Feldes und 
einer Abschatzung der Orientierung des Generators (7), und 
IS 2) Venivenden der Abschatzungen der Entfemung aus Schritt 1 ) und des bekannten Ortes des Fekisensors 

(1) zum Berechnen einer Abschatzung des Ortes des Generators (7); 

3) Venvenden des geschatzten Ortes des Erzeugers (7) aus Schritt 2) und der Messungen des Feldes bei 
den Sensoren (1) zum Berechnen einer neuen Abschatzung der Orientierung des Erzeugers (7); 

4) iteratives VViederholen der Schritte 1) bis 3), wobei in Schritt 1) zum Verbessern der Abschatzungen des 
20 Ortes und der Orientierung des Erzeugers (7) die neue Abschatzung der Orientierung des Erzeugers von dem 

vorausgehenden Schritt 3) venwendet wird. 

8. Vorrichtung zum Bestimmen des Ortes und/oder der Orientierung eines Feldsensors (7). mit: 

2S mehreren Fetderzeugern (1 ), wobei jeder mehrere zusammengestellte felderzeugende Elemente (3) enthalt, 

einer Einschattein richtung (5), die fur jeden Felderzeuger (1) so betrieben werden kann, daB sie mehrere 
felderzeugende Elemente (3) zum Erzeugen eines Feides, innertnalb dessen der Ort und die Orlentlemng des 

Sensors (7) bestimmt werden konnen, einschalten kann, 

einer an den Feldsensor (7) angeschlossenen MeQeinrichtung (8) zum Messen und Ausgeben eines Para- 

30 meters des Feldes bei dem Sensor (7), und 

einer Steuerelnrichtung (6) zum Steuern der Einschalteinrichtung (5), wobei die Steuereinrichtung (6) eine 
Speichereinrichtung zum Speichem von durch die MeBeinrichtung (6) ausgegebenen Parametern und eine 
Recheneinrichtung zum Berechnen des Ortes und der Orientienjng des Sensors (7) aus den gespelcherlen 
Messungen des Feldes enthalten, 

35 wobei die Vorrichtung dadurch gekennzeichnet ist, da3 die Recheneinrichtung folgendes enthalt: 

eineOrtsabschatzungseinrichtung, die so betrieben werden kann, daB sie aus den gespeicherten Messungen 
der Felder und einer Abschatzung der Orientierung des Sensors (7) eine Abschatzung der Entfemung des 
Sensors (7) von jedem Felderzeuger (1) berechnen kann und unter Venwendung der Abschatzungen der Ent- 
femung und des bekannten Ortes des Felderzeugers (1) eine AbschStzung des Ortes des Sensors (7) be- 

40 rechnen kann; 

eineOrientierungsabschatzungseinrichtung, die so betrieben werden kann, daB sie aus der Abschatzung des 
Sensorortes durch die Ortsabschatzungsein richtung und die gespeicherten Messungen des Feldes eine neue 
Abschatzung der Orientierung des Sensors (7) berechnen kann; . 

eine Einrichtung zum iterativen Vertsessem der Abschatzung des Ortes und der Orientierung des Sensors (7). 
4S wobei dtese Einrichtung beim Gebrauch wiederholt sequentiell die Oris- und Orientierungsabschatzungsein- 

richtung verwendet, bei dem die AusgangsgroBen der Orientierungsabschatzungseinrichtung bei jeder Itera- 
tk>n von der OrtsabschStzungseinrichtung als die neue AbschStzung der Sensororientierung verwendet wer- 
den. 

so 9, Vorrichtung zum Bestimmen des Ortes eines Feldsensors (7) nach Anspruch 8. bei dem der Feldsensor (7) ein 
einzelnes Feldsensorelement enthalt und bei dem der durch das einzelne Feldsensorelement abgetastete Para- 
meter des Feides von der Orientierung des Feldsensorelements innerhalb des FeMes abhangt. 

10. Vorrichtung zum Bestimmen des Ortes eines Feldsensors (7) nach Anspmch 8 oder Anspruch 9 mit drel Felder- 
ss zeugern (1), von denen jeder drei zueinander orthogonale zusammengestellte felderzeugende Elemente (3) ent- 
halt 

11. Vorrichtung zum Bestimmen des Ortes eines Objekts relativ zu einem Felderzeuger (1), wobei die Vorrichtung 
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enthalt 

einen Felderzeuger (1) mit mehreren felderzeugenden Elementen (3), 

eine Einschalteinrichtung (5) zum sequentiellen Einschaften jedes der felderzeugenden Elemente (3) zum 
s Erzeugen eines Feides, innerhaib dessen das Objekt lokalisiert werden kann, 

einen an dem Objekt anzubringenden Sensor (7) mit mehreren zusammengestellten Feldsensoreiementen 
(A, B, C), die jeweils einen Parameter des Feides abtasten konnen, 

eine Me3einrichtung (6) zum Messen und Ausgeben des Parameters des durch jedes Feldsensorelement (A, 
B. C) abgetasteten Feides. wenn jedes felderzeugende Element (3) eingeschaltet 1st, und 
10 eine Steuereinrichtung (6) zum' Steuem der Einschalteinrichtung (5), zum Spelchem der durch die MeBefn- 

richtung (8) ausgegebenen Parameter und zum Berechnen diss Ortes des Sensors (7) relative zu dem Felder- 
zeuger (1), 

wobel die Vorrichtung dadurch gekennzetchnet ist. da3 die Steuereinrichtung (6) enthalt: 

IS 

eine Gewlchtungseinrichtung zum Gewichten jedes der von dem Sensor (7) abgetasteten Parameter durch 
einen richtungsabhangigen Gewichtungsfaktor fOr jedes felderzeugende Element (3), so da3 die berechnete 
Feldstarke B gleich der Feldstarke B ist, die bei dem Sensor (7) existieren wQrde, falls die Achse des Feides 
zum Sensor (7) gerichtet ware, 
^ eine Einrichtung zum iterativen Andem der richtungsabhangigen Gewichtungsfaktoren, bis ein ^4aximatwert 

der Feldstarke get unden ist, und zum Besttmmen der Richtung des Sensors (7) von dem Felderzeuger (1 ) bis 
zu einem gewunschten Genauigkeitsgrad, und 

eine Recheneinrichtung zum Berechnen der Entfernungdes Sensors (7) von dem Erzeuger (1) und daherdes 
Ortes des Sensors (7) relativ zu dem Erzeuger (1). ; 

2S 

12. Vorrichtung zum Besttmmen des Ortes eines Objekts nach Anspruch 11 , bei dem der Felderzeuger (1 ) drei Zuein- 
ander orthogonale felderzeugende Elemente (3) enthSIt und der Sensor (7) drei zueinander orthogonale Feldsen- 
sorelemente (A, B, C) enthalt. 

30 13. Vorrichtung nach irgendeinem der Anspruche 6 bis 12, die so beschaffen ist, d.a3 sie h der Chirurgie, Therapie 
Oder Diagnose eines menschlk:hen Oder tierischen Korpers venA^endet wird. 

14. Endoskopiesystem zur Verwendung in der Endoskopie eines menschltchen Oder tierischen Korpers, wobel das 
Endoskopiesystem eine in irgendeinem der Anspruche 8 bis 12 beanspruchte Vorrichtung enthalt. . 

35 

15. Vorrichtung zum Uberwachen des Ortes eines Endoskops innerhaib eines menschlichen Oder tierischen Korpers, 
wobel die Vorrichtung versehen ist mit einem Fekjsensor (7), 

mehreren Felderzeuger (1), wovon jeder mehrere felderzeugende Elemente (3) enthalt, 
40 einer Einschalteinrichtung (5), die fur jeden Felderzeuger (1) so betrieben werden kann, daQ sie jedes der 

mehreren felderzeugenden Elemente (3) zum Erzeugen eines Feides. innerhaib dessen der Ort des Endo- 
skops bestimmt werden kann, einschalten kann, 

einer an den Feldsensor (7) angeschlossenen Me8einrichtung (8) zum Messen und Ausgeben eines Para- 
meters des Feides bei dem Sensor (7), und 
45 einer Steuereinrichtung (6) zym Steuem der Einschalteinrichtung (5), zum Speichern der durch die Meflein- 

richtung (8) ausgegebenen Parameter und zum Berechnen des Ortes des Sensors (7) relativ zu den mehreren 
Felderzeugem (1), wobei die Vorrichtung dadurch gekennzeichnet ist, daB die Steuereinrichtung (6) so be- 
schaffen Ist. daB sie den Ort des Sensors (7) relativ zu den mehreren Felderzeugem (1 ) mit dem in irgendeinem 
der AnsprOche 1 bis 7 beanspruchten Verfahren berechnet. 

50 

16. Vorrichtung zum Bestimmen des Ortes und/oder der Orientierung eines Sensors (7) nach irgendeinem der An- 
spruche 8. 9 Oder 10, bei dem einer Oder mehrere der Sensoren (7) an Teilen eines Handschuhs befestigt sind. 
wobei der Ausgang der Vorrbhtung Daten bereitsteltt. die den Ort und die Orientierung der jeweiligen Teile der 
Hand eines Tragers des Handschuhs definieren. 

55 

17. Vorrichtung nach Anspruch 16, bei dem einer Oder mehrere Sensoren (7) an Fingern des Handschuhs befestigt 
sind. 
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Revendlcations 

1. Proc6d6 de determination de la localisatton et/ou de Torientation d'un capteur de champ (7) par rapport k une 
pluralit6 de g6n6rateurs de champ (1) de localisation connue, chaque g6n6rateur de champ (1) comprenant une 
5 plurality rf6l6ments de g6n6ratlon de champ situ6s au mfime emplacement (3), chaque 6l6ment de g6n6ration 

(3) dtant activ6 et une mesure 6tant faite au niveau du capteur de champ (7) du champ respectif g6n6r6 par 
I'6l6ment de g6n6ratlon (3). le proc§d6 6tant caract^risd par les 6tapes consistant ^ : 

1) pour chaque g6n6rateur de champ (1). calculer ^ partir des mesures du champ g6n6r6 par chacun de ses 
10 dldments de gdn^raticn (3) et d'une estimation de rorientation du capteur (7), une estimation de la dlstan6e 

ehtre ce g^ndrateur de champ particulier (1) et le capteur (7), et 

2) utiliser les estimations de la distance provenant de P6tape 1), et la localisation connue des g6n6rateMrs do 
champ (1 ) pour calculer una estimation de la localisation du capteur (7). 

3) utiliser la localisation estim6e du capteur (7) provenant de P6tape 2) et les mesures du champ au niveau 
IS du capteur (7) pour calculer una nouvelle estimation de i'orientation du capteur (7), 

4) r^pdter les stapes 1) ^ 3) de fagon iterative. I'dtape 1) utiiisant la nouvelle estimation de I'orientation du 
capteur provenant de r^tape pr6c6dente 3), afin d'amdiiorer les estimations de la localisation et de i'orientation 
du capteur (7). 

20 2. Procddfi selon la revendlcation 1 , dans lequel chaque 6l6ment de g6n6fation (3) de chaque g6n6rateur de chanap 

(1) est actlvd en sequence de sorts qu'un seul element de generation (3) est active k un instant donne. 

3. Precede selon la revendication 1, dans lequel au moins deux elements de generation (3) sont activ6s simultane- 
ment, chacun de ces au moins deux elements de generation (3) etant active k une frequence differente de rautre 

2s ou des autres elements de generation actives simultanement (3). 

4. Precede selon la revendication 1 . dans lequel k l'6tape (3), la localisation estimee du capteur (7) calcuiee k retape 

(2) et les mesures du champ au niveau du capteur (7) sont utilisees pour calculer, pour chaque generateur de 
champ (1 ), la direction du champ respectif au niveau du capteur (7) et, k partir de ces directions de champ, calculer 

30 une nouvelle estimation de I'orientation du capteur (7). 

5. Precede selon la revendication 1 ou la revendication 4, dans lequel k I'etape (1 ). I'estimation inltiale de I'orientation 
du capteur (7) est prise comme etant I'orientation qui a ete calcuiee par I'application Immedlatement pr6cedente 
du proc6de ou, s'il s'agit de la premiere application du precede, Taxe du capteur (7) est suppose initlalement 6tre 

35 oriente en direction du generateur de champ respectif (1). 

6. Precede de determination de la localisation d'un capteur de champ (7). comprenant une pluralitd d*616ment8 de 
detection de champ situ6s au mfime emplacement (A, B, C) par rapport k un g6n6rateur de champ (1 ) comprenant 
une pluralite d'eiements de generation de champ silues au mdme emplacement (3). chaque element de generation 

40 de champ (3) etant active tour k tour et la valeur de I'intensite du champ resultant etant mesuree au niveau du 

capteur de champ (7), lintensite du champ mesuree etant d6pendante de la localisation du capteur (7) k I'interieur 
du champ respectif, le precede etant caracterise par les etapes consistant k : 

1 ) calculer, en utiiisant toutes les valeurs d'intensite de champ mesur6es et une estimation de la direction du 
4S capteur (7) par rapport au generateur de champ (1 ), un facteur de pond6ration dependant de la direction pour 

chaque element de generation de champ (3) de sorte que I'intensite de champ calcuiee B soit 6gale k I'intensite 
du champ B qui existerait au niveau du capteur (7) si I'axe du champ etait oriente en direction du capteur (7), 

2) modifier de fa^on iterative les facteurs de pond6ration dependant de la direction afin de maximiser B et de 
determiner ainsi avec un niveau desire de precision la direction du capteur (7) par rapport au generateur de 

so champ (l),et 

3) utiliser les valeurs mesurees de I'intensite du champ pour calculer la distance du capteur (7) par rapport 
au generateur de champ (1 ) et de ce fait, k partir de la direction du capteur (7) de retape 2), la localisation du 
capteur (7) par rapport au gendrateur de champ (1). 

ss 7. Precede de detemnination de la localisation et/ou de I'orientation d'un generateur-de champ (7) par rapport k une 
pluralite de capteurs de champ (1 ) de localisation connue, chaque capteur de champ (1 ) comprenant une pluralite 
d'eiements de capteur de champ situ6s au mdme emplacement (3), le generateur (7) etant active et une mesure 
etant faite au niveau des elements de capteur de chanrp (3) du champ respectif genere par le generateur (7), le 
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proc6dd 6tant caractdris6 par les Stapes consistant k : 

1 ) pour chaque capteur de champ (1 ), calculer a partir des mesures du champ g6n6r6 au niveau de chacun 
de ses 616ments de capteur (3) et d'une estimation de I'orientation du gSnSrateur (7), une estimation de la 

5 distance entre ce capteur partlculier(1 ) et le g6n6rateur (7). et 

2) utiliser les estimations de distance provenant de PStape 1 ). et la localisation connue des capteurs de champ 
(1) pour calculer une estimation de la localisation du gSnSrateur (7). 

3) utiliser la localisation estlmde du g§nSrateur (7) provenant de f^tape 2) et les mesures du champ au niveau 
des capteurs (1) pour calculer une nouvelle estinnation de i'orientation du gSndrateur (7). 

10 4) rSpdter les Stapes 1) ^ 3) de fa9on iterative, PStape 1) utilisant la nouvelle estimation de I'orientation du 

gSnSrateur provenant de I'Stape prScddente 3). afin d'amSliorer les estimations de la localisation et de rorien- 
tation du gSnSrateur (7). 

8. Oispositif destinS d determiner la localisation et/ou I'orientation d'un capteur de champ (7) comprenant : 

15 

une plurality de gSnSrateurs de champ (1) comprenant chacun une pluraiitS d'SISments de gSnSration de 
champ situSs au mdme emplacement (3), 

un moyen d'activation (5) pouvant 6tre mis en oeuvre pour chaque g6n6rateur de champ (1) afin d'activer 
chacun de la plurality des elements de gSndratbn de champ (3) de mani^re k produire un champ ^ I'intSrieur 
20 duquel la localisation et I'orientation du capteur (7) peuvent §tre ddterminSes, 

un moyen de mesure (8). reli6 au capteur de champ (7), afin de mesurer et de foumir en sortie un param6tre 
du champ au niveau du capteur (7), et 

un moyen de commando (6) afin de commander le moyen d'activation (5), le moyen de commando (6) com- 
prenant un moyen de memorisation destjn6^ m6moriser des paramfetres fournis en sortie par le moyen de 
ss mesure (8). et un moyen de calcul destine k calculer, d partir des mesures memorisSes des champs, la loca- 

lisation et Torientation du capteur (7), le dispositif etant caracterisd en ce que le moyen de calcul comprend : 

un moyen d'estimation de localisation qui peut §tre mis en oeuvre pour calculer, k partir des mesures 
m6moris6es des champs et d'une estimation de I'orientation du capteur (7). une estimation de la distance 
30 du capteur (7) par rapport ^ chacun des gdnSrateurs de champ (1 ) et, en utilisant les estimations de 

distance et la localisation connue des g^nSrateurs de champ (1 ). calculer une estimation de la localisation 
du capteur (7), 

un moyen d'estimation d'orientation qui peut dtre mis en oeuvre pour calculer, k partir de Testimation de 
la localisation du capteur par le moyen d'estimation de localisation et des mesures m6moris6es du champ, 

3S une nouvelle estimation de I'orientation du capteur (7), 

un moyen destine d ameiiorerde fa^on iterative Testimation de la localisation et de I'orientation du capteur 
(7), lequel moyen, en utilisation, utilise de fa^on repetitive les moyens d'estimation de localisation et 
d'orientation en sequence, dans lequel la sortie du moyen d'estimation d'orientation est adoptee par le 
moyen d'estimation de localisation en tant que nouvelle estimation de I'orientation du capteur h chaque 

40 iteration. 

9. Dispositif destine d determiner la localisation d*un capteur de champ (7) sebn la revendicatlon 8, dans lequel le 
capteur de champ (7) comprend un 616ment de detection de champ unique et le paramMre du champ detecte par 
ledit element de detection de champ unique depend de I'orientation de I'diement de detection de champ k I'lnterieur 

4S du champ. 

10. Dispositif destine k determiner la localisation d'un capteurde champ (7) selon la revendicatlon 8 ou la revendicatlon 
9, comportant trois genSrateurs de champ (1), chacun comprenant trois 6l6ments de generation de champ situ6s 
au mdme emplacement mutuellemenl orthogonaux (3). 

50 

11. Dispositif destine k determiner la localisation d'un objet par rapport k un g6nerateur de champ (1). le dispositif 
comprenant 

un g6n6rateur de champ (1) comprenant une pluralite d'6iements de generation de champ (3), 
55 un moyen d'activation (5) destine k activer s6quentiellement chacun des elements de generation de champ 

(3) afin de produire un champ k I'lnterieur duquel Tobjet peut dtre localise, 

un capteur (7), destine h 6tre fixe h I'objet. comprenant une pluralite d'eiements de detection de champ situes 
au m§me emplacement (A, B. C), chacun pouvant detector un paramfetre du champ, 
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un moyen de mesure (8) destind k mesurer et k fournir en sortie le param^tre du champ ddtectd par chaque 
^Idment de dStectbn de champ (A. B. C) iorsque chaque 6l6ment de gdndration de champ (3) est activ6. et 
un moyen de ccmmande (6) destind k commander le moyen d'activation (5), k m^moriser les paramMres 
foumis en sortie par le moyen de mesure (8) et k calculer la localisation du capteur (7)..par rapport au g^n§* 
5 rateurde champ (1), 

le dispositif 6tant caract^risd en ce que le moyen de commando (6) comprend 

un moyen de pond^ration destind k ponderer chacun des paramdtres, ddtect^s par le capteur (7). par un 
facteur de pond6ration dependant de la direction pour chaque 6l6ment de gdniration de champ (3), de sorte 
que I'intensitd caiculde du champ B soit 6gale k Tintensitd du champ B qui existerait au niveau du capteur (7) 
10 si i'axe du champ dtait orients en direction du capteur (7), 

un moyen destind k modifier detain iterative les facteurs de ponddration dependants de la direction jusqu*^ 
ce qu'une valeur nnaximum de Tintensitd du champ soit trouv^e. et k determiner la direction du capteur (7) par 
rapport au g^nSrateur de champ (1) avec un niveau d^sir^ de precision, et 

un moyen de calcul destine k calculer la distance du capteur (7) par rapport au gen^rateur (1 ) et de ce fait la 
IS localisation du capteur (7) par rapport au g^nerateur (1 ). 

12. Dispositif destine k determiner la localisation d'un objet seton la revendication 11, dans lequet le generateur de 
champ (1 ) comprend trois elements de generatbn de champ mutuellement orthogonaux (3) et le capteur (7) com- 
prend trois elements de detection de champ mutuellement orthogonaux (A, B, C). 

so 

13. Dispositif selon Tune quelconque des revendlcations 8^12, agence pour une utilisation en chiairgie, theraple ou 
diagnostic d'un corps humain ou animal. 

14. Systems d'endoscopie destine k etre utilise lors de I'endoscopie d'un corps humain ou animal, le systems d'en- 
25 doscopie comprenant un dispositif selon Tune quelconque des revendications 8 £i 12. 

15. Dispositif destine ^ contrdter la position d'un endoscope k i'interieur d'un corps hunnain ou animal, le dispositif 
comprenant un capteur de champ (7). 

30 une pluralite de generateurs de champ (1) comprenant chacun une pluralite d'eiements de generation de 

champ (3), 

un moyen d'activation (5) pcuvant etre mis en oeuvre pour chaque generateur de champ (1) afin d'activer 
chacun de la pluralite des elements de generation de champ (3) de maniere k produire un champ k I'interieur 
duquel la localisation de I'endoscope peut 6tre determinee, 
35 un moyen de mesure (8) relie au capteur de champ (7), afin de mesurer et de fournir en sortie un pararhetre 

du champ au niveau du capteur (7), et 

un moyen de commando (6) destine k commander le moyen d'activation (5), memoriser les parametres foumis 
en sortie par le- moyen de mesure (8). et calculer la localisation du capteur (7) par rapport k la pluralite des 
generateurs de champ (1), le dispositif etant caractdrise en ce que le moyen de commando (6) est agence 
40 pour calculer la localisation du capteur (7) par rapport k la pluralite des generateurs de champ (1) grkce au 

precede selon Tune quelconque des revendications 1^7. 

16. Dispositif destine k determiner la localisation et/ou forientation d'un capteur (7) selon Tune quelconque des re- 
vendications 6, 9 ou 10. dans lequel un ou plusieurs desdits capteurs (7) sont fixes k des parties d'un gant, la 

45 sortie dudit dispositif procurant des donnees qui definissent la localisation et Porientation des diverses parties de 
la main du porteur du gant 

17. Dispositif selon la revendication 16, dans lequel un ou plusieurs capteurs (7) sont fixes aux doigts du gant. 

so 
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Fig. 5. 
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Fig. 7. 
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Fig. 8. 




31 



EP 0 655 138 B1 

Fig. 10. 
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Fig . 12 
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Fig . 13 
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